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ABSTRACT 

We present here, for the first time, a 2D study of the overshoot convective mecha- 
nism in nova outbursts for a wide range of possible compositions of the layer underlying 
the accreted envelope. Previous surveys studied this mechanism only for solar com- 
position matter accreted on top of carbon oxygen (C-O) white dwarfs. Since, during 
the runaway, mixing with carbon enhances the hydrogen burning rates dramatically, 
one should question whether significant enrichment of the eject a is possible also for 
other underlying compositions (He, O, Ne, Mg), predicted by stellar evolution models. 
We simulated several non-carbon cases and found significant amounts of those under- 
lying materials in the ejected hydrogen layer. Despite large differences in rates, time 
scales and energetics between the cases, our results show that the convective dredge 
up mechanism predicts significant enrichment in all our non-carbon cases, including 
helium enrichment in recurrent novae. The results are consistent with observations. 
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1 INTRODUCTION 

Almost all classical and recurrent novae for which reliable 
abundance determinations exist show enrichment (relative 
to solar composition) in heavy elements and/or helium. It 
is now widely accepted that the source for such enrichment 
is dredge up of matter from the underlying white dwarf to 
the accreted envelope. A few mechanisms for such mixing 
were proposed to explain the observations : Mixing by a dif- 
fusion layer, for which diffusion during t he accretion phase 
builds a layer of mixed ab u ndanc es (iPrialnik &, Kovetzl 



dl994h: 



Shankar et al. 



ing t he runaway itself | Wooslev 
(1 19921) ; IShankar fe Arnett 

dl995h: [Glasner et al. I (Il99<t): iKercek et al. I (11998 



Glasner &; Livne 
1999h 



iGlasner et al. I (j2007T ) : [Casanova et al. I (|2Qlol . I20TT 



dl985h : llben et al. I (|l99ll . 

iFuiimoto fe Iben I (|l992j V)~ Mixing by shear insta- 



iKovetz &; Prialnik 



bility in duced by d ifferen ti al rotation during the ac c retion 
phase (|Durisen Il977l; [K ippenh ahn & Thomas! Q .978; 
MacDonald Ill983l : iLivio fe Truran Ill987l : iKutter fe Sparks I 
19871. 119891 ): Mixing by shear gravity waves breaking 
on the white dwarf surface in which a resonant inter- 
action between large-scale shear flows in the accreted 
envelope and gravity waves in the white dwarf's core 
can induce m i xing o f heavy elements into t he en velope 
(|Rosner et al. I l200ll : lAlexakis et~al~l 12001 l2004h and 
finally - Mixing by overshoot of the convective flow dur- 



In this work we focus on the last of these mechanisms 
that proved efficient for C-0 white dwarfs. Mixing of carbon 
from the underlying layer significantly enhances the hydro- 
gen burning rate. The enhanced burning rate drives higher 
conve ct ive fluxes, inducing more m i xing I Glasner &; Livne I 
(1 1995T ): [Glasner et al. I (|l997l . I2007T ) : ICasanova et al. I (|2010l . 
l2011a lb). Therefore, the fact that the underlying layer is 
rich in 12 C is a critical feature of all the overshoot convec- 
tive models that have been analyzed up to this work. Ac- 
cording to the theory of stellar evolution for single stars, we 
expect the composition of the underlying white dwarf to be 
C-0 for masses in the ran ge 0.5 - 1.1 MQ and ONe(Mg) fo r 
more massive whit e dwar fs[Gil-Pons & Garcia Berro ( 2Q0ll): 
iDomingues et al. I <|2002h : iGarcia-Berro et al. I (|2002l ). Ob- 
ser vat ions show enric hment in helium, CNO, Ne, Mg 
and heavier eleme nts ( Starrfield, Sparks & Truran I [l986; 
iGehrz et al~l Il998l . 120081 : llliadis et al. I l2002h . For recur- 
rent novae, helium enrichment can achieve levels of 40 — 
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50% (IWebbink et al. Ill987l : lAnupama and Dewanganll2000l : 
iDiaz et al. Il2010h . High helium abundances can simply be 
explai ned as the ashes of hydrogen burning during the run- 
away ([Hernanz.M. &; Jose. J. Il2008l ) , but one can not exclude 
the possibility that the source of He enrichment is dredge up 
from an underlying helium layer. We therefore found it es- 
sential to study nova outbursts for which the composition of 
the underlying layer is different from C-O. The models stud- 
ied here extend the work we publish in the past. As a first 
step we study here the runaway of the accreted hydrogen 
layer on top of a single white dwarf, changing only its com- 
position. Having a fixed mass and radius we can compare the 
timescales, convective flow, energetics and dredge up in the 
different cases. A more comprehensive study which varies 
the white dwarf's mass with compositions is left to future 
work (CO, ONe(Mg) or He rich). 

The present study is limited to 2D axially symmet- 
ric configurations. The well known differences between 2D 
and 3D unstable flows can yield uncertainties of few per- 
cents on our results, but can not change the general trends, 
as previous studies showed reasonable agreement between 
2D and 3D simulations with regard to integral quantities, 
although larger differen ces persist in the local structure ( 
ICasanovaet a l~l (fioTTbl ^. We therefore regard our present 
results as a good starting point to more elaborated 3D sim- 
ulations. In the next section we describe the tools and the 
setup of the simulations. In subsequent sections we describe 
the results for each initial composition and then summarize 
our conclusions. 



2 TOOLS AND INITIAL CONFIGURATIONS 

All the 2D simulations presented in this work start from a ID 
hydrostatic configurations, consisting of a 1.147M© CO core 
with an outer layer of 10 -4 Mq composed of CO, ONe(mg) 
or helium, according to the studied case as we explained in 
the introduction. The original core was built as a hydro- 
static CO politrop that cooled by evolution to the desired 
central temperature (2 x 10 7 K). The ID model evolves using 
Lagrangian coordinates and does not include any prescrip- 
tion for mixing at the bottom of the envelope. The core ac- 
cretes matter with solar abundance and the accreted matter 
is compressed and heated. Once the maximal temperature 
at the base of the accreted envelope reaches a temperature 
of 9 x 10 7 K, the whole accreted envelope 3.4 x 10~ 5 M© and 
most of the underlying zone 4.7 x 10 -5 M© , are mapped 
onto a 2D grid, and the simulations continue to r unaway 
and b eyond using the 2D hydro code VULCAN-2D iLivne I 
(|l993h . This total mass of 8.1 x 10~ 5 M© is refered as the 
total computed envelope mass. Using same radial zoning in 
the ID grid and in its 2D counterpart, the models preserve 
hydrostatic equilibrium, accurate to better than one part 
in ten thousand. Since the configurations are unstable, non 
radial velocity components develop very quickly from the 
small round-off errors of the code, without introducing any 
artificial initial perturbation. 

Further computational details of the 2D simulations 
are as follows. The inner boundary is fixed, with assumed 
zero inward luminosity. The outer boundary follows the ex- 
pansion of the envelope taking advantage of the arbitrary 
Lagrangian- Eulerian (ALE) semi- Lagrangian option of the 



Table 1. Parameters of the Simulated Initial Configurations 



Model 


U nderlying 


Tmax 


Qmax 


remarks 


ml2 




2.05 


4.0 


ID 


ml2ad 


CO 


2.45 


1000.0 




ml2ag 


He 








ml2dg 


He 






T base = 1.5 X 10 8 K 


ml2al 


O 








ml2bl 


O 






T base = 1.05 x 10 8 K 


ml2cl 


O 






T base = 1.125 x 10 8 K 


ml2dl 


O 


2.15 


82.4 


T base = 1.22 x 10 8 K 


ml2kk 


O 






24 Mg rates 


ml2jj 


O 


2.46 


1000.0 


24 Mg rates+ T base = 1.125 x 10 8 X 



Tmax'i maximal achieved temperature [10 8 Kelvin]. Qmax', maxi- 
mal achieved energy generation rate [10 42 erg/sec] 

VULCAN code whereas the burning regions of the grid at 
the base of the hydrogen rich envelope are purely Eulerian. 
More details are presented in Glasner et al. (2005, 2007). The 
flexibility of the ALE grid enables us to model the burning 
zones at the bottom of the hydrogen rich envelope with very 
delicate zones in spite of the post-runaway radial expansion 
of the outer layers. The typical computational cell at the 
base of the envelope, where most of the burning takes place, 
is a rectangular cell with dimensions of about lAkm x 1.4/cra 
. Reflecting boundary conditions are imposed at the lateral 
boundaries of the grid. Gravity is taken into account as a 
point source with the mass of the core, and the self gravity 
of the envelope is ignored. The reaction network includes 
15 elements essential for the hydrogen burning in the CNO 
cycle, it includes the isotopes: 1 H, 3 He, 4 He 7 Be, 8 B, 12 C, 

13 C 13 N 14 N 15 N 14q 15q 16q 17q and 17 R 



3 RESULTS 

In order to compare with ID models and with previous stud- 
ies we present here five basic configurations: 

1) The outburst of the ID original model without any 
overshoot mixing. 

2) An up to date model with an underlying C-0 layer. 

3) A model with a Helium underlying layer. 

4) A model with an underlying O(Ne) layer. 

5) A toy model with underlying 24 Mg. 

This model demonstrates the effects of possible mixing 
of hydrogen with 24 Mg on the runaway. (In a realistic model 
the a mounts of 24 Mg in the core are exp e cted to be a fe w per- 
cents (|Gil-Pons fc Garcia Berrol (|200lh : ISiess I <|2006h ), but 
highe r amounts can be found in th e very outer layers of the 
core (|Garcia-Berro fc Iben 1 (|l994h ). 

The computed models are listed in Table [TJ In the next 
sections we present the energetics and mixing results for each 
of the models in the present survey. 



3.1 The ID original model 

The initial (original) model, composed of a 1.147M© de- 
generate core, accretes hydrogen rich envelope at the rate of 
1.0 x 10 -10 Mq /year. When the temperature at the base of 
the accreted envelope reaches 9 x 10 7 K the accreted mass 
is 3.4 x 10 _5 M©. We define this time as t=0, exceptional 
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test cases are commented in Table U Convection sets in for 
the original model a few days earlier, at t = — 10 6 sec, when 
the base temperature is 3 x 10 7 K. The ID convective model 
assumes no overshoot mixing; therefore, convection has an 
effect only on the heat transport and abundances within 
the convective zone. Without the overshoot mixing burning 
rates are not enhanced by overshoot convective mixing of 
CO rich matter. As a reference for the 2D simulations, we 
evolved the ID model all through the runaway phase. The 
time to reach the maximal energy production rate is about 
2400sec, the maximal achieved temperature is 2.05 x 10 8 K, 
the maximal achieved burning rate is 4.0 x 10 42 erg /sec and 
the total nuclear energy generated up to the maximum pro- 
duction rate (integrated over time) is 0.77 x 10 46 erg. 

3.2 The C-O underlying layer 

We summarize here the mai n results of the unde rlying C- 
O case (computed already in Idasner et al. I (|2007h and re- 
peated here), which is the most energetic case. A comparison 
of the history o f the burning rate (Fig. [1]) with Figure 3 in 
Idasner et al. I (|2007h confirms that our current numerical 
results agree with those of our earlier publication. In this 
figure (Fig. [TJ we also present the amount of mixing at vari- 
ous stages of the runaway. The main effects of the convective 
underlying dredge up are: 

• The convective cells are small at early stages with mod- 
erate velocities of a few times 10 6 cm/sec. As the energy 
generation rate increases during the runaway, the convec- 
tive cells merge and become almost circular. The size of the 
cells is comparable to the height of the entire envelope, i.e. 
a few pressure scale heights. The velocity magnitude within 
these cells, when the burning approaches the peak of the 
runaway, is a few times 10 7 cm/sec. 

• The shear convective flow is followed by efficient mixing 
of C-0 matter from the core to the accreted solar abundant 
envelope. The amount of C-0 enrichment increases as the 
burning becomes more violent and the total amount of mix- 
ing is above 30% (Fig.HJ. 

• Mixing induces an enhanced burning rate, relative to 
the non mixing ID case, by more than an order of magni- 
tude. The maximum rate grows from 4.0 x 10 42 erg /sec to 
1000.0 x 10 42 erg /sec. The enhanced rates raise the burn- 
ing temperature and shorten the time required to reach the 
maximal burning rate. The maximal achieved temperature 
increases from 2.05 x 10 8 K to 2.45 x 10 8 K and the rise time 
to maximum burning decreases from 2440sec to 140sec. The 
total energy production rates of the ID and the 2D simula- 
tions are given in Fig. [1] The enhanced burning rate in the 
2D case will give rise, at later stages of the outburst, to an 
increase in the kinetic energy of the eject a. Unfortunately, 
since the hydro solver time step is restricted by the Courant 
condition we can not run the 2D models through to the 
coasting phase. A consideration of the integrated released 
energy at the moment of maximal burning rate reveals that 
the burning energy grows from 0.77 x 10 46 erg in the ID 
model to 1.45 x 10 46 erg in the 2D model, a factor of 2. 

Another interesting feature of the 2D C-0 simulations 
is the appearance of fluctuations, observed during the initial 
stages of the runaway (Fig. [J). Such fluctuations are not ob- 
served in the ID model. The fluctuations are a consequence 



8 44 




(a) 




Mass/Tot_mass 



(b) 



— 1D no mixing 

— 2D CO layer 




Log(Time) [sec] 

(c) 

Figure 1. The CO model a - the logarithm of the energy gen- 
eration rate as function of time, b- the total abundance of the 
CNO elements as function of the fractional mass coordinate of 
the computed envelope (see section 2) at various times that are 
indicated on the plot, c - the log of the total energy generation 
rate as function of the logarithm of the time coordinate for the 
ID model and for the 2D CO model. 



of the mixing of the hot burning envelope matter with the 
cold underlying white dwarf matter. The mixing has two ef- 
fects. The first is cooling as we mix hot matter with cold 
matter. The second is heating by the enhancement of the 
reaction rate. It is apparent that in this case, after a small 
transient, the heating by enhanced reaction rates becomes 
dominant and the runaway takes place on a short timescale. 
For other underlying compositions the effect is a bit more 
complicated. We discuss this issue in the next subsections. 




Time [sec] 

Figure 2. Log of the total energy production rate for all the 
helium models compared to the ID model. The model with initial 
temperature higher than the default temperature of 9 x 10 7 K was 
shifted in time by 1760.0 seconds. 



3.3 The Helium underlying layer 

The helium enrichment, observed in recurrent nova (without 
enrichment by heavier isotopes) and in other classical nova 
was mentioned in the past as an obstacle to the un derlying 
convective overshoot mechanism (|Livio &; Truranl (| 19871 )). 
Helium is the most abundant end product of the hydrogen 
burning in nova outbursts and it does not enhance the hy- 
drogen burning in any way. In recurrent novae, helium may 
be accumulated upon the surface of the white dwarfs. If so, 
can the dredge-up mechanism lead to the observed helium 
enrichment ? We examine this question using the case ml2ag 
in Table [U Energetically, as expected, the model follows ex- 
actly the ID model (Fig. [2]). The slow rise of the burning 
rate in this case makes the 2D simulation too expensive. 
To overcome this problem we artificially 'jump in time' by 
jumping to another helium model at a later stage of the 
runaway, m which the ID temperature is 1.22 x 10 8 K. The 
rise time of this model is much shorter, and by adjusting its 
time axis to that of the ID model the two curves in Fig. [2] 
may be seen to coincide. The fluctuations of the 2D curve 
are absent in its ID counterpart, both because the latter has 
no convection and because the ID simulation is performed 
using implicit algorithm with much larger time steps. The 
precise way in which the 2D evolution agrees with the ID 
evolution increases our confidence in the validity of the 2D 
simulations. 

The convective flow is indeed moderate but overshoot 
mixing is observed at a certain level (Fig. [3J. In the first 
(slower) phase (ml2ag in Table [T]) the level of mixing is 
small, converging to about 10%. In the later (faster) phase 
(ml2dg in Table [1]) the mixing rate increases in step with 
the increasing burning rate and the higher velocities in the 
convective cells. However, as the rates are still low relative 
to the C-0 case, the added amount of mixing is again about 
10%, summing up to a total amount of about 20%. Since we 
begin with matter of solar composition, the total He mass 
fraction at the end of the second phase exceeds 40%. 

The color maps of the absolute value of the velocity in 
the 2D models at different times along the development of 




0.2 h 
0.1 j- 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
Mass/Tot_mass 

(b) 

Figure 3. Abundance of the helium as function of the fractional 
mass coordinate of the computed envelope (see section 2) at var- 
ious times that are indicated on the plot, a- the model with the 
default base temperature of 9 x 10 7 X, b - the model with base 
temperature of 1.22 x 10 8 X (times shifted by 1760 seconds, the 
2D time is in brackets). 

the runaway are presented in Fig. 2] In the first slower phase 
(model ml2ag Table [1]), the burning rate is low and grows 
mildly with time. The convective velocities are converging 
to a value of a few 10 6 cm/sec and the cell size is only a 
bit larger than a scale height. In the second (faster) phase 
(model ml2dg Table[T|), the burning rate is somewhat higher 
and is seen to grow with time. The convective velocities are 
increasing with time up to a value of about 1.5 x 10 7 cm/sec. 
The convective cells in the radial direction converge with 
time to an extended structure of a few scale heights. 

3.4 The ONe(Mg) underlying layer 

The rate of proton capture by oxygen is much slower and 
less energetic than the capture by carbon but it still has an 
enhancing effect relative to the ID model without mixing 
(Fig. [5}. This can be well overstood once we notice that for 
the initial temperature of the 2D model i.e. 9.0 x lO 7 ^ the 
energy generation rate by proton capture by oxygen is more 
than three orders of magnitude lower than the energy gen- 
eration rate by carbon capture (Fig. [6}. We can also observe 
that the energy generation rate for capture by oxygen stays 
much smaller than the energy generation rate for capture by 
carbon for the entire range of temperature relevant to nova 
outbursts. For this range, we also notice that the capture 
by 20 Ne is lower by an order of magnitude than the capture 
rate by 16 O. Being interested only in the energetics, convec- 
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(e) 

Figure 4. Color map of the convective speed for the underlying helium layer. The times are: model ml2ag (a)- t=50 sec,(b)- t=200 sec, 
(c)- t=300 sec; model ml2dg (d)- t=40 sec, (e)- t=140 sec , (f)- t=230 sec 



tive flow and mixing by dredge up, we choose to seperate 
varibles and study the case of a pure 16 O underlying layer 
as a test case (ignoring any possible 24 Mg that is predicted 
by evolutionary codes). 

We computed the model for an integrated time of 300 
seconds (about one million time steps). The trend is very 
clear and we can extrapolate and predict a runaway a bit 
earlier and more energetic than the ID case. Again, we could 



not continue that 2D simulation further due to the very low 
burning rates and a small hydrodynamical time step. As be- 
fore, we computed three different phases, where each of them 
starts from a different ID model along its evolution. The 
maximal ID temperatures at the base of the burning shell, 
when mapped to the 2D grid, are : 1.05 x 10 8 K, 1.125 x 10 8 K 
and 1.22 X 10 8 K respectively (Tabled . All three phases start 
with a transient related to the buildup of the convective flow. 
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Time [sec] 

Figure 5. Log of the total energy production rate for all the 
oxygen models compared to the ID model and to the CO model. 
The models with initial temperature higher than the default tem- 
perature of 9 x 10 7 K were shifted in time in order to produce a 
smooth continuous line. 




Temperature [K] 



Figure 6. Log of the energy generation rate for proton capture 
on: 12 C, 16 , 20 Ne, and 24 Mg, the rate is calculated for p = 
1000.0 gr/cc. 



In Fig. [5] we shifted the curves of burning rates in time in a 
way that permits a continuous line to be drawn. Along this 
continuous line, evolution proceeds faster towards a runaway 
than the ID burning rate, also shown in Fig. [5] 

We computed the last phase, with an initial base tem- 
perature of 1.22 x 10 8 i^ , for 350 sec until it reached a maxi- 
mum. The maximal achieved temperature is 2.15 x lO 8 ^ and 
the maximal achieved burning rate is 82.4 x 10 42 erg /sec. As 
expected, this case lies somewhere between the ID case and 
the C-0 model. The convective flow resembles the C-0 case, 
a significant feature being the strong correlation between 
the burning rate and the convective velocities. Most impor- 
tantly for the case of an underlying oxygen layer, dredge up 
of substantial amounts of matter from the core into the en- 
velope occurs in all our simulations. The trends are about 
the same as for the C-0 models (Fig. [7}. The correlation of 
the amount of mixing with the intensity of the burning is 
easily observable. 
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Figure 7. Abundance of the CNO elements as function of the 
fractional mass coordinate of the computed envelope (see section 
2) at various times that are indicated on the plot, a - the model 
with the default base temperature of 9 x 10 7 K, b - the model 
with base temperature of 1.22 x 10 8 K. 



3.5 Approximate models for Mg underlying 
layer 

Nova outbursts on massive ONe(Mg) white 
dwarfs are expected to be e nergetic fast nova 
dStarrneld. Sparks Truran I Il98& iGehrz et al.l Il998l . 
2008: llliadis et al. 1 120021 ). The problem we face is whether, 
in absence of the enhancement by 12 C, the overshoot mixing 
mechanism can generate such energetic outbursts by mixing 
the solar abundance accreted matter with the underlying 
ONe(Mg) core. Based upon examination of the energy 
generation rate of proton capture reactions (p, 7) on 12 C, 
16 0, 20 Ne , and 24 Mg the results shown in Fig. ED make it 
evident that only 24 Mg can compete with 12 C in the range 
of temperatures relevant to nova runaways. Therefore, in 
spite of the fact that the abundan ce of 24 Mg in the core 
sums u p to o n ly a few percent ([Gil-Pons &; Garcia Berrol 
([200 lh : ISiessl (|2QQ6h ). the high capture rate might com- 
pensate for the low abundances and play an important 
role in the runaway. Furthermore, previous studies show 
that in the outer parts of the core, the parts important for 
our study, 24 M g is more abundant a n d can represent up 
to about 25% (Garcia-Berro & Iben ] (|l994h ). Restricting 
ourselves to the reaction network that includes only 15 
elements we assume, as a demonstration, an artificial case 
of a homogeneous underlying layer with only one isotope 
( 24 Mg). For this homogeneous layer model we replaced 
the energy generation rate of proton capture by 16 O with 




— Modified rates - Mg mixture - used for 2D model 

— Exact rates - Mg mixture 

— Exact rates - CO mixture 



le+08 1.5e+08 2e+08 2.5e+08 

Temperature [Kelvin] 

Figure 8. Log of the total energy production rate for p = 1000.0 
gr/cc of a mixture that contains 90% solar matter mixed with 
10% of 24 Mg. Red: 15 elements net used for the 2D model; Blue: 
The rates given by a full net of 216 elements. The Black line gives 
the rates of 90% solar matter mixed with 10% of CO core matter 
(see text). 
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Figure 9. Log of the total energy production rate for all the 
magnesium models compared to the CO model the ID model 
and the oxygen model. The model with initial temperature higher 
than the default temperature of 9 x 10 7 K is not shifted in time 
(see text). 



the values of proton capture by 24 Mg Fig. [6] To check our 
simplified network, we present in Fig. [8] the rates computed 
by a 216 elements network and our modified 15 elements 
rates, both for a mixture of 90% solar matter with 10% 
24 Mg. The difference is much smaller than the difference 
inside the big network between this mixture and a mixture 
of 90% solar matter with 10% C-0 core matter. Therefore 
our simplified network is a good approximation regarding 
energy production rates. 

The crossing of the curves in Fig.[8]reveals a striking and 
very important feature - at temperatures less than 1.3 x 10 8 K 
a mixture of 10% carbon and 90% solar compostion burns 
roughly ten times faster than a mixture of 10% magnesium 
with the same solar composition gas. Above that tempera- 
ture the rates exchange places and magnesium enhancement 
dominates C-0 enhancement. This emphasizes the impor- 
tance of a proper treatment of the effects of the 24 Mg abun- 
dance in explosive burning on ONe(Mg) white dwarfs. In a 
future work, we intend to study more realistic models with 
the inclusion of a detailed reaction network. 

In Fig. [9] we present the total burning rates in our toy 
model, together with the rates of previous models. 

As can be expected, the enhancement of the burning in 
this toy model, relative to the ID model, is indeed observed. 
However, the development of the runaway although faster 
than in the underlying 16 O model is still much slower than 
the rise time of the C-0 model. This result is easily overstood 
via the discussion above, as the initial burning temperature 
in the model is only ( 9x 10 7 K). At that temperature, the en- 
ergy generation rate for proton capture on 24 Mg is lower by 
almost three orders of magnitude relative to the energy gen- 
eration rate by proton capture on 12 C. The rates are about 
the same when the temperature is 1.3 x 10 8 K and from there 
on the magnesium capture rate is much higher. In order to 
demonstrate that this is indeed the case we calculated an- 
other 2D magnesium model in which the initial maximal ID 
temperature was 1.125 xl0 8 K. The rise time of this model is 
very short, even shorter than the rise time of the C-0 model. 
One should regard those two simulations as two phases of 



one process - slow and fast. The maximal achieved temper- 
ature is 2.45 x 10 8 K and the maximal achieved burning rate 
is 1000.0 x 10 42 erg /sec, similar to the C-0 case. 

In order to better understand the convective flow for 
the case with underlying 24 Mg we generated color maps of 
the absolute value of the velocity (speed) in the 2D mod- 
els at different times along the development of the runaway 
(Fig. Eg. The two 24 Mg cases show extremely different be- 
havior. In the first case (model ml2kk Table [1]), the burning 
rate is low and it grows mildly with time. The convective 
velocities are converging to a value of a few 10 6 cm/sec and 
the cell size is only a bit bigger than a scale height. In the 
second case (model ml2jj Table [1]), the burning rate is high 
and it grows rapidly with time. The convective velocities 
are increasing with time up to a value of a few 10 7 cm/sec. 
The convective cells in the radial direction converge to a 
structure of a few scale heights. 

In accordance with our previous cases, the magnesium 
toy model dredge up substantial amounts of matter from the 
core to the envelope. There is a one to one correlation with 
the convective velocities. The amount of mixing at the slow 
initial stages (model ml2kk) are small and tend to converge 
to a few percents. The amount of mixing at the late fast 
stages grows rapidly with time (Fig. [TT]) . We present here 
only the general trend detailed results will be presented in 
a forthcoming study. 



4 CONCLUSIONS 

We present here, for the first time, detailed 2D modeling 
of nova eruptions for a range of possible compositions be- 
neath the accreted hydrogen layer. The main conclusion to 
be drawn from this study is that significant enrichment 
(around 30%) of the ejected layer, by the convec- 
tive drege-up mechanism, is a common feature of 
the entire set of models, regardless of the composi- 
tion of the accreting white dwarf . On the other hand, 
the burning rates and therefore the time scales of the run- 
away depend strongly on the composition of the underlying 



8 Ami Glasner et al. 




3e+D7 
. 2,7e+07 
. 2,4e407 

1.8e407 
l,5e+07 



(a) 





3e+D7 R 



(c) 
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Figure 10. Color map of the convective speed for the underlying magnesium layer. The times are: model ml2kk (a)- t=40 sec, (b)- 
t=60 sec, (c)- t=70 sec; model ml2jj (d)- t=20 sec, (e)- t=30 sec ,(f)- t=40 sec 



layers. There is also a one to one correlation between the 
burning rate, the velocities in the convective flow, and the 
amount of temporal mixing. Therefore, second order differ- 
ences in the final enrichment are expected to depend on the 
underlying composition. Specific results for each case are as 
follows : 

a) Since the energy generation rate for the capture of 
protons by 12 C is high for the entire temperature range pre- 



vailing both in the ignition of the runaway and during the 
runaway itself, the underlying carbon layer accelerates the 
ignition and gives rise to C-0 enrichment in the range of the 
observed amounts. 

b) For the densities and temperatures prevailing in nova 
outbursts helium is an inert isotope. Therefore, it does not 
play any role in the enhancement of the runaway. Neverthe- 
less, we demonstrate that once the bottom of the envelope 
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Figure 11. Mixing for the underlying magnesium 
layer, abundance of the CNO elements as function of the 
fractional mass coordinate of the computed envelope (see section 
2) at various times that are indicated on the plot, a - the model 
with the default base temperature of 9 x 10 7 K, b - the model 
with base temperature of 1.125 x 10 8 K. 



a source for burning enhancement by overshoot mixing is 
24 Mg. From our demonstrating toy model, we can speculate 
that even small amounts of 24 Mg present at the high stages 
of the runaway can substantially enhance the burning rate, 
leading to a faster runaway with a significant amount of 
mixing. The relationship between the amount of 24 Mg in 
the ONe(Mg) core, the steepness of the runaway, and the 
amount of mixing in this case are left to future studies. 
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is convective, the shear flow induces substantial amounts 
of mixing with the underlying helium. The eruption in 
those cases is milder, with a lower burning rate. For recur- 
rent nova, where the timescales are too short for the diffu- 
sion process to play a significant role, the observed helium 
enrichment favor the underlying convection mechanism as 
the dominant mixing mechanism. Future work dealing with 
more realistic core masses (1.35-1.4 solar masses) for recur- 
rent novae will give better quantitative predictions that will 
enable us to confront our results with observational data. 

c) The energy generation rate for the capture of a proton 
by 16 O is much lower than that of the capture by 12 C for 
the entire temperature range prevailing in the ignition of the 
runaway and during the runaway itself. Underlying oxygen, 
whenever it is present, is thus expected to make only a minor 
contribution to the enhancement of the runaway. As a result 
the time scale of the runaway in this case is much larger 
than that of the C-0 case. Still, the final enrichment of the 
ejecta is above 40%, (Fig. [7J). The energy generation rate 
by the capture of a proton by 20 Ne is even lower than that 
of capture rate by 16 O. We therefore expect 20 Ne to make 
again only a minor contribution to the enhancement of the 
runaway, but with substantial mixing. 

d) Nova outbursts on massive ONe(Mg) white dwarfs 
are expected to be energetic fast nova. In this survey we 
show that for the range of temperatures relevant for the nova 
runaway the only isotope that can compete with the 12 C as 
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